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Abstract 

This article concerns the time-dependent Hartree-Fock (TDHF) approximation of 
single-particle dynamics in systems of interacting fermions. We find that the TDHF 
approximation is accurate when there are sufficiently many particles and the initial many- 
particle state is any Gibbs equilibrium state for noninteracting fermions (with Slater de- 
terminants as a special example). Assuming a bounded two-particle interaction, we obtain 
a bound on the error of the TDHF approximation, valid for short times. We further show 
that the error of the TDHF approximation vanishes at all times in the mean field limit. 



1 Introduction 

The time-dependent Hartree-Fock (TDHF) equation is a nonlinear Schrodinger equation de- 
signed to approximate the evolution of an n-electron system. The TDHF equation was first 
written down by Dirac, both as a system of n coupled Schrodinger equations for occupied or- 
bitals, and as an integro-differential equation for the "density matrix," i.e., the integral kernel 
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in atomic units, where V^xt is the external potential energy and V{r) = 1/r is the Coulomb 
interaction potential. The the last term on the right-hand side of (Q) is the "exchange" term. 

The Coulomb potential, however, is not amenable to the techniques of this article because it is 
unbounded (this case being dealt with in In this article, we consider interaction potentials 
given by a bounded function V{\x — y\), or, more generally, any bounded, symmetric, two-body 
operator V (not necessarily a multiplication operator). The one-particle energy operator will 
be denoted by L, the interaction energy operator for a single pair of particles will be denoted by 
V, and the total energy operator for a system of particles will be the sum of all single-particle 
energies and all pair energies. Although the number of particles does not change under the 
dynamics just described, we prefer to formulate the dynamics on a fermion Fock space so that 
we may consider initial states of indeterminate particle number. We are going to show that 
quasifree initial states enhance the accuracy of the TDHF approximation. 

Let L be a self-adjoint operator on a Hilbert space H, and let y be a bounded Hermitian 
operator on EI(8>IHI that commutes with the transposition operator U defined by U{x®y) = y®x. 
We are going to discuss the dynamics whose Hamiltonian H on the fermion Fock space Fh is 
written as 



in second quantized form. We will analyze the solutions of the von Neumann equation 



which is the evolution equation for the density operator on Fock space in the Schrodinger 
picture of quantum dynamics (in units of time and energy for which ^ = 1). 

We will see that © leads to the following equation for the single-particle (number) density 
operator Niit): 
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where [V, Af2{t)];i denotes the partial trace of the two-particle operator [V, A/2(t)]. Equation 
dD) for A/i(t) is not "closed" since its right hand side involves the two-particle density operator 
A/2(i)- The TDHF approximation to A/i(t) is the solution of the initial value problem 




[L,F{t)] + [V, F{tf'2A2].. 

Mp(o)) 
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where A2 is the orthogonal projector of EI EI onto the subspace of antisymmetric vectors. The 
existence and uniqueness of solutions of © were established in [71 for the case where V is a 
bounded operator, and in [TU[ |H] for the case where V is a. Coulombic interaction. 

The TDHF equation (jSJ is obtained by closing the single-particle equation (jH) with the 
Ansatz 



at all times. The relation (jH)) holds for pure states corresponding to Slater determinants, and 
also for Gibbs densities. However, even supposing that A/2(0) satisfies the interaction V 
is likely to introduce "correlations" in N'2{t), that is, departures from (jU)), and ignoring those 
correlations in the TDHF equation requires justification. 

We are going to prove that the absence of correlations is self-perpetuating in the mean field 
limit. Theorem 16.21 states that if M2{0) satisfies (jH)) then Af2(t) asymptotically satisfies © as 
the number of particles n tends to infinity and the interaction strength is scaled as 1/n. In 
this scaling, the force exerted on each individual particle by the n — 1 other particles is of 0(1) 
as n — ^ cxo. This was called the "mean-field scaling" by H. Spohn in his fundamental review 
paper where he derives the time-dependent Hartree equation. In the appendix we prove 
an important special case of Theorem 5.7 of |14j as a corollary of our Theorem 16.21 

Technically, we rely on the trace norm approach used in to derive the Hartree equation. 
We first published our derivation of the TDHF equation in the mean field limit in |Hj for initial 
states of fixed particle number such as Slater determinants. The main advance of this article 
is that the initial states need not be Slater determinants; the TDHF approximation should 
work equally well (or badly) for any quasifree initial state. Also, in this article we are not 
only interested in the mean field limit, and we derive the error bound of Theorem 16. II for the 
unsealed problem. 

Throughout this paper, n denotes a fixed number of particles, denotes the the number 
operator on the Fock space over EI, and A/i denotes the number density operator on EI. 

2 Definitions and notation 

Consider a quantum particle whose Hilbert space is EI, i.e., a particle which, in isolation, would 
constitute a system whose (pure) quantum states are represented by the rank-one orthogonal 



ATs = {Afi^Afi)2A2 
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projectors on some Hilbert space H. The set of quantum states available to a system of n 
particles of this kind is determined by their "statistics," i.e., whether the particles are fermions, 
bosons, or distinguishable. If the particles are fermions, a (pure) state of a system of n of them 
is represented by a rank-one projectors onto vectors in the antisymmetric subspace H'^") of the 
tensor power space H®"^. If the number of fermions in the system is not fixed, the appropriate 
Hilbert space is the fermion Fock space 

Fh = h(°) © e(^) © e(2) © h(3) © • • • . (7) 

The possibility of a zero-particle state is accommodated by M^^\ a one-dimensional space 
spanned by the vacuum vector. We denote the number operator on Fh by N. 

Let n„ denote the group of permutations of {1,2,..., n}. For each n G n„, a unitary operator 
on H®" may be defined by extending 

Ut,{xi ^ X2 ■ ■ ■ Xn) = ® a;^-i(2) ® ■ ■ ■ ® x^-^n) 

to all of H^". The operator 

= ^ 5^ sgn(7r)[/^ (8) 
7ren„ 

is the orthogonal projector with range H("). If X\^X2-, ■ ■ ■ iXn is an orthonormal system in a 
Hilbert space H, then the vector 

^\ An{Xx®X2® ■ ■ ■ ®Xn), (9) 

is a unit vector in EI'^"\ called a Slater determinant. 

We will consider many-particle states that are represented by density operators on Fh. We 
will only consider density operators that commute with and such that N"^D is trace class 
for all m G N. For such densities D, one can define the reduced density operators Nm{D) of all 
orders m: 

If T is a trace class operator on EI^"\ and m < we will use the subscript notation 
to denote the partial trace of T of order m, a trace class operator on H*^™'-'. (This operator is 
defined unambiguously thanks to the symmetry of T considered as an operator on H®*^.) If 
T is a density operator on EI*^"\ the operator T-m is known as the m-particle reduced density 
operator JH] and it is used to determine the expected values of the m-particle observables. Let 
Z) be a density operator on Fjj that commutes with A^. Then 

OD 

D = (10) 

n=0 

where each Dn is a nonnegative trace class operator on M^'^\ Assuming that 

oo 

Tr(Ar™L)) = ^n'"Tr(D„) < oo, (11) 

n=0 
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we may define the m order reduced density operator 

oo , 

Af^D) = J2l '-^D^.m, (12) 

^ — ^ Ira — ml! 

n=m ^ ' 

where Dnvm = {Dn):m is the partial trace of Dn- The reduced density Mm serves to describe 
the m-particle correlations in a system of many particles (see Section 6.3.3 of 



3 Dynamics and the BBGKY hierarchy 

We now define the dynamics Q and the reduced dynamics 

Let if be a self-adjoint operator on H, and let be a bounded Hermitian operator on EI (8> EI 
that commutes with the transposition operator U{x®y) = y®x. For 1 < j < n, let Lj denote 
the operator 

j—l times n—j times 

/ (S) ■ ■ ■ ® / ®L® / (g) • ■ ■ (g) / 

on H®" (the value of n > j is not explicit in the notation Lj but it will always be clear from 
context). For 1 < i < j < ra, let ?7(jj) denote the permutation operator on H®" that transposes 
the i*^^ and j^^ factors of any simple tensor xi (g ■ ■ ■ x„, and let 

= f/(l.)f/(2,)(V^®/^"-')f/(2,)f/(l.) 

(again, the domain H*®" of Vij will always be clear from context). For each n, define the 
operators 

n 

l<i<j<?i 

on H^'^) (these operators are defined on all of H®'^ but we are only considering their restrictions to 
the invariant subspace H^")). The Hamiltonian operator which we had formally represented 
above by Q, is the direct sum H = ^ Hi"^^ defined on the domain 

V{H) = |x = ©a;„eFH: ^ ||ii(")x„||^ < oo}. 

n 

This operator is closed and self-adjoint (see Section 6.3.1 of [9,) and —iH is the generator of 
the strongly continuous group 

oo 
n=l 
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of unitary operators on Fh, where W} = exp ( — itH'^^'^Y The djTiamics corresponding to Q 
are given by the group 

WtiD) = WtDW.t (13) 

of isometries of the space of Hermitian trace class operators on Fe- (See Proposition 3.4 of 
[7j for a proof that groups of isometries defined in this way are strongly continuous.) It is 
convenient to have some notation for the free part of the dynamics, so we define 

W,(")(T) = [/.("^Tf/if (14) 

with U^"'^ = exp {-itL^'^^). 

The dynamical equation for the m*'^ order reduced density Mm can be derived from ()13p if 
the density D satisfies the moment condition (lllj). The details of the derivation are provided 
in gj. 

Proposition 3.1 Let Ut and he as defined in and M'J\) . 

Suppose that D is a density operator on Fe of the form D = ©-D„ such that ill]) holds for 
somemEN. Let Afmif) denote Afm(yVt{D)) . Then Afm(t) satisfies 




Mn+l(s)].„ ds . 



These equations for the reduced density operators are known as the BBGKY hierarchy. The 
first equation of the hierarchy is equation (g} in integral form. 



4 The TDHF hierarchy 

The existence and uniqueness of mild solutions of the TDHF equation ^ is established in [Zj. 
There it is shown that the integral equation 

F{t) = W«F(0) - ^ fut\ [V, F{sr'2A2\.^ds (15) 

Jo 

has a unique solution F{t) for any Hermitian trace class operator -F(O). Define J-'i{t) = F{t) 
and, for m > 1, define 

J-„(t) = Fitf^mlAm. (16) 

We proceed to derive equations for the J-'mif) from (fT^. 
First, set G{t) = U^}}F{t), so that 

G{t) = F(0) - t fu^}] [V, F{s)^^2A2].^ds. 
Jo 



Now apply the product rule (in integral form) to G{t) 



m „t 
,=1 ^0 



F(0) 



-1 ^7/(1) 



m 



t m 
(m) 



r r 

I / Wi^5^ krn+1, F(.)^-+l (/-[/(,,„+!)) 

^0 L 



(is. 



Apply W^*-™''' to both sides of the preceding equation to obtain 



ds. 



Multiply both sides of the last equation by m\Am on the left, noting that lA'f^\x)Ar, 
Us^\XAm) and that commutes with ^ Vj^m+i- 



J^^t) = W^)^^(O) - t / F(s)«"+i(/-%„+i))m!A, 



ds. 



U(m,m+i))TTT'^-Am, the last equatiou may be rewritten 



Since (m + l)!v4m+i = (/ - ^{i,m+i) " • ■ 

„( .11. 

We call these equations for the J-'m{t) the TDHF hierarchy. 



ds 



(17) 



5 Estimates 

In this section we have collected some estimates used in the proofs of Theorems 16 . II and 16 . 21 The 
first two of the following propositions are stated here without proof, for the facts are known, 
and the reader may find proofs of them in We use the notation || ■ || for the operator norm 
and II ■ 111 for the trace norm. 



Proposition 5.1 If D is a density operator on Fe that commutes with N and such that 
Tr{ND) < oo, then the operator norm of N'i{D) is not greater than 1. 
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Proposition 5.2 If T is a Hermitian trace class operator then 

||T®"n!A„||^ < ||t||". (18) 

Proposition 5.3 The trace norm of the last term in ^17}) does not exceed 

2m(m- 1)||V|| ||F(0)|| ||i^(0)||7t . 

Proof: The trace norm of the last term in (|17|) is bounded by 

2m{m-l) {Vm-i,m+iU(m,m+i){J^m{s) ^ F{s))},^ ds (19) 
Jo ■ 1 

thanks to the symmetry of J^m- It can be verified that 
whence (fTUj) is bounded by 

2m(m-l) f \\V\\ \\F{s)\\ ||j^„(s) 
Jo 

Now ||JSn(s)||i < ||i^(s)||^" by Proposition O since Tm = F®"^m\A^. Furthermore, ||F(s)||i = 
||F(0)||i and ||F(s)|| = ||F(0)|| for all s > by Proposition 4.3 of ^. Thus we arrive at the 
bound stated in the Proposition. □ 

6 Accuracy of the TDHF approximation 

In this section we will compare the single-particle density operator Mi{t) to its approximation 
by the solution of the TDHF equation 

F{t) = ui'^'FiO) - I y"*wS[^, F{sr^2A,],^^ds 

F(0) = M(0). (20) 

We can control the distance between A/i(t) and F{t) in trace norm when the initial many- 
particle state is gauge-invariant quasifree. Theorems 16 . II and 16 . 21 show how the accuracy of the 
TDHF approximation is enhanced when the initial state involves many uncorrelated fermions. 

We consider the class of initial states that have density operators D of the form (jiup such 
that Tr(A^L') < oo and 

Mn{D) = M(Z})®"n!A„ (21) 
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for all n G N. These relations characterize the gauge-invariant quasifree states of the CAR 
algebra having a trace class single-particle reduced density operator Afi P . 

There are two important examples of such initial states: Slater densities and Gibbs grand 
canonical equilibrium states. "Slater densities" are those of the form 

n—l times 

D = 0©---©0 ©o©--- , 

where is an ra-particle Slater determinant Gibbs equilibrium states are obtained as 

follows: let H be the single-particle Hamiltonian and dr{H) its functorial extension to the 
fermion Fock space viewed as the exterior algebra of H. In other words, 

n>0 

where dVn{H) is the restriction to H^"^ of the operator 

" J — 1 times n—j times 

/©...© / (gH® /©...©/ 

J=l 

Whenever /5 > is such that e~^^ is trace-class on H, the Gibbs equilibrium state at inverse 
temperature /3 with chemical potential fi is defined by the density operator proportional to 
exp{—(3dr{H — yu/)) (see Proposition 5.2.22 of [HI for more details on Gibbs states). 

Theorem 6.1 Let D be the density operator on Fjj of a gauge-invariant quasifree state with 
finite expected particle number, i.e., with Tr{ND) < oo. Let Afiit) denote Ni{Wt{D)) , where 
Wt is the dynamics with two-particle interactions defined in 

Let F{t) be the solution of the TDHF equation \2U\] . 

Let T denote (211^/1111^111)-^ Then 

\\m)-F{t)l < (22) 

for t < T. 

The proof of this theorem is postponed until the end of this Section. 

Unfortunately, the bound ()22|) on the error of the TDHF approximation is valid only when 
t < T and we have no explicit bounds for larger t. In principle, the estimate of Theorem 16.11 
could be used to establish that any effect observed in a TDHF simulation before the critical time 
r reflects a true effect of the interaction. However, in the numerical tests we have conducted 
so far, very little appears to happen before the critical time r, and we fear that the estimate of 
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Theorem 16.11 might not prove generally useful. Nonetheless, Theorem 16.11 does show that the 
error of the TDHF approximation is less than one might expect when the initial condition is 
an uncorrelated many-fermion state, for even at short times t < r the left-hand side of ()22|) is 
proportional to ||F(0)||i = ||A/i(-D(0))||i prima facie, not bounded independently of ||A/i||i. 

The improvement in accuracy of the TDHF approximation when the initial state is uncor- 
related is even more evident in the "mean field scaling." It is in this spirit that we are about 
to introduce a coupling constant A into the interaction term of the many-particle system and 
consider a scalings where A times the average particle number tends to or remains bounded. 
We do not discuss the physical significance of such scalings; we consider them only so that we 
may more easily express how the accuracy of the TDHF equation is affected by uncorrelated 
initial data. 

For each value of the parameter A > 0, consider the Hamiltonian 

Hx = ^^(j|L|z)ajaj + A V"|zj)a|,a|ajaj. (23) 

i,j i,j,k,l 

If the initial density operator represents a gauge-invariant quasifree state with finite expected 
particle number, then Proposition 13.11 implies that all reduced number density operators exist 
and satisfy 

/•t 

•^^ l<i<j<m j=l 

AfxmiO) = AfmiDxiO)). (24) 



The TDHF equation corresponding to (j23j) is 

F,{t) = W«F,(0) - \^ [V, F,isr'2A,],^ds 

Jo 

Fx{0) = AfxiiO). (25) 

Theorem 6.2 Let {Dx}^^^ be a family of density operators on Fe that represent gauge-invariant 
quasifree states with finite expected particle number. Let Nxmif) be the solution of with 
initial condition A^Am(O) = Nm{Dx) and let Fx{t) be the solution of the TDHF equation \25]) . 
Let J'xmit) = Fxitr^^mlA^. 

1) // lim A||A/'ai(0)||i = then, for each fixed t > and m eN, 



x^o 



\^^Xrr.{t)-J'xn.{t)y\\^fxl\\T = 0(A). 
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2) // limsup A||7Vai(0)||i < oo then 

\im\\^Ut) - J'Ut)\\J\\^x4T = 

A — ^0 

for allt>0 and all m G N. 

Remark on the persistence of the interaction in the mean field hmit 

Although we have assumed only that the interaction ^ is a bounded and symmetric two- 
particle operator, we are mainly interested in the case where it is a bounded multiplication 
operator. In this case one may bound the exchange term in the TDHF equation as follows. 
The exchange term in (Q) is the sum of two products: the product of the integral operators with 
kernels V{\x — z\)F{x,z) and F{z,y), and the product of the integral operators with kernels 
F{x,z) and V{\y — z\)F{z,y). The trace norm of each of these products is bounded by the 
product of the Hilbert-Schmidt norms of its factors. Since the operator norm of a fermionic 
single-particle operator is less than or equal to 1 (viz. Proposition 15. ip 

\\F\\hs = VTr(F2) < ^/WMf\\ < VWh 

and it follows that the trace norm of the exchange term is bounded by 2||F||i||y||j;^oc. Thus the 
contribution of the exchange term is not much larger than the error of the TDHF approximation 
itself! To see this more clearly, let us return to Statement 1 of Theorem 16.21 Squeezing the 
most we can out of its proof informs us that 

\\^fM{t) - F,{t)l = 0{X'\\F,\\l) 

if limA||FA||i = 0. (It is due to the fact that the j = term of the second series on the 

right-hand side of the bound in Lemma f6. II below vanishes when m = 1.) But the part of Fx{t) 
due to the exchange term is 0(A||Fx||i) by the above estimate, and A||Fa||i already tends to 
zero. 

Knowing that the exchange effect vanishes in the mean field limit, one might wonder whether 
Theorem 16.21 is trivial. The theorem states that the error between the true single-particle 
density and its TDHF approximation tends to zero under certain conditions — but perhaps 
this so simply because the effect of the interaction disappears in the limits we have taken? At 
least the "direct" part of the interaction does not disappear in the limit, for it survives as the 
nonlinear term in the time-dependent Hartree equation as shown in the following corollary of 
Theorem 16.21 

Corollary 6.1 Let ip e L^(M^) have norm one, and for n e N let \E'„(a;,t) G L^(]R^"') be the 
solution of the Schrddinger equation 

d 1 " 1 

i — <i/n{x,t) = --J2^x,'^n{x,t) + Xjl)"^ n{x,t) 

j=l i<j 

^„(x,0) = V^(xi)7/'(x2)---V^(a;„) (26) 
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with V{x) bounded. Let pn{x,y,t) denote the integral kernel of the single-particle density oper- 
ator, normalized to have trace 1 rather than n: 



Pn{x,y,t) = / ■■■/ '^ri{y,Z2,...,Z.n,,t)'^n{x,Z2,...,Zn,t)dZ2---dZn. (27) 

Let p^{x,y,t) denote the solution of the time- dependent Hartree equation 

d f 

'i^P*{x,y,t) = -UA.^- Ay)p^{x,y,t) + / [V{\x - z\) - V{\y - z\)]p^{z, z,t)dz p^{x,y,t) 

Ot .hi3 



p4x,y,0) = ij{y)ij{x) . (28) 

Then Pn{t) converges in trace norm to p^:{t) at each fixed t > as n — > oo. 

The proof of this corollary is given in the appendix. 

Proof of Theorems 16.11 and 16.21 

Let N'm{t) be as in Proposition and let Tm{t) satisfy the TDHF hierarchy. In the 
hypotheses Theorems 16.11 and 16.21 we suppose that F{fS) = J\fi{D{0)), but for now let us only 
assume that 

||F(0)||<1 and l|F(0)lK = l|M(/^(0))||i. (29) 

The trace norm of A/i(t) is independent of t, and we shall denote it simply by ||A/i||i. Assuming 
(EHjl . the bound of Proposition 15.31 is itself bounded by 



mini 



1)2\\V\\ llMll^t. (30) 



Subtracting equations (fTTjl from the BBGKY equations of Proposition 13. II and using (jHUjl leads 
to the estimates 



\^/mit) ~ ^m{t)\\, < |K™(0)- J-„(0)||^ + m(m-l)||y||(^2||M|r + l|A^mllij^ 

+ m2||V|| / \\j\fm+i{s) - T.m+i{s)\\^ds . 
Jo 



'0 

Iterating this estimate k times, one obtains 

k / . ■ ,\ k 



+ '^"t TTT / / ■■■/ \\-^m+k+i{s) - J^m+k+i{s)\\-^dsdtk---dti 

(31) 
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with C = 2\\V\\ and 



hm. = m 



Ar^(O)- J-„(0)||^ 

(m-l)||V||(2||M||r+||A/'„||i). (32) 



To make use of these estimates we need some control over the size of the integrand in . We 
will assume that 

IIAT^lli < \m\T (33) 

for all m, for this bound holds for gauge-invariant quasifree state with finite expected particle 
number by (PT|) and Proposition 15.21 Note that (jHHj) holds independently of time, since the 
dynamics conserve particle number. With the bound (j33|) . the last term on the right hand side 
of may be bounded by 



which tends to as A; tends to infinity if m fixed and C||A/i||it < 1. Furthermore, assuming 
(j33p ■ we can bound bm of ^ by lC\\Afi\\ ^m{m — 1) and establish the following lemma: 

Lemma 6.1 Suppose that D is a density operator on Fe of the form D = ®Dn, such that 
holds for all m G N. Let Wt he as defined in / f73|) and let Mmif) denote A/'m(Wt(Z^)). Let F{t) 
he the solution of a TDHF equation / f73)) whose initial condition -F(O) satisfies /f^). and let 
Tmif) he as in 17^) . Then, with C = 2||V"||, 



\\^m{t)-' 


^ m (^) 1 


ll 




|M| 


1 m 
ll 



j=0 

oo 



when C\Mx\\t < 1. 



In the hypotheses of Theorems 16.11 and 16.21 the initial data for the exact dynamics are 
assumed to be gauge-invariant quasifree states D (or Dx) with finite expected particle number, 
and the initial data for the corresponding TDHF equations are assumed to be F{0) = A/i(-D) 
(or Fx{0) = J\fi{Dx)). Thus, the requirements and in Lemma IHTD are satisfied under 
the hypotheses of Theorems 16.11 and 16.21 

Theorem 16 . II follows from Lemma [6.1 1 since ^^^(0) = J^m(O) for all m by ^IT^ . 

To prove Theorem 16. 21 we apply Lemma f6. II to the many-particle system (|T^ and the TDHF 
equation fT5|) with W in place of V. 
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Statement 1 of the theorem comes easily: Mxm{0) = J-'xmi^) for all m, and Lemma IHTD implies 
that 

with C = 2\\V\\ when ACIIA/ai Iht < 1. Statement 1 follows since lim AIIA/aiIIi = 0. 

The proof of Statement 2 requires the fuller version of the inequality in Lemma l6.1[ Since 
()33|) and (j^^ are satisfied at any time s > if they are satisfied initially, Lemma 16.11 implies 
that 



\MUs + At) - rx,a(s + At)\\^ ~ \\{Mx),„+j{s)-{J',U+,(s)lfm + j-l\ 



Al 1 



for any s > as long as At < (AC||A/ai||i)^^. Since u = Climsup A||A/ai||i is finite by 

A-»0 

hypothesis, taking the lim sup of both sides of (jH^ shows that 

lim \\AfxUt) - :FxUt)\\ J WxiWT = VmGN (35) 



holds at time t = s + At if it holds at time t = s and At < 1/u. Since f)35|) holds at t = 0, an 
inductive argument proves that it holds at all times t > 0. □ 



7 Appendix: the time-dependent Hartree equation 

Here we prove Corollarv 16.11 

The initial condition of ()26|) is not available to n fermions in M^, for it is in extreme violation 
of the Pauli Exclusion Principle. To derive (j^Hj) from our theorem about fermions we will 
introduce an auxilliary space £^(N) to allow Pauh exclusion to hold while the spatial part of 
the n-particle wavefunction is permitted have the form 'il!{xi)ip{x2) ■ ■ -ip^Xn)- We are going to 
apply Statement 2 of Theorem 16 . 21 where the single-particle Hilbert space EI = L'^(R^) ® 

Let ei, 62, . . . be an orthonormal sequence in £^(N) and let Di/„(0) be the orthogonal projector 
on Fe whose range is the span of the n-particle Slater determinant 

An {{tp ® ei) ® (^Z- ® 62) ® ■ ■ • ® (^A ® en)) . 
It is helpful to rearrange factors and write 

Di/n{0) = RniP^^n Sn)Rn , 
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where P^-gm denotes projection onto the span of t/'®"', and 5"^ denotes the rank-one orthogonal 
projector on (?{H)®'^ whose image is the span of the Slater determinant formed from ci, . . . , e„, 
and Rn is the unitary transformation from L'^{M?)®'' ® ^{n)®'' to {L'^{M?) that 
rearranges the factors of any simple tensor thus: 

^n((/l ® f2® ■ ■ ■ ® fn) ® {Si® S2® ■ ■ ■ ® S„)) = (/i ® Si) (g) ■ ■ ■ ® (/„ ® S„) . 

If ^nif) is the solution of then P^!^{t) — the projector whose range is the span of '^nif) — 
satisfies the von Neumann equation 

d 1 " 1 

j=l i<j 

wherein Myij denotes the multiplication operator Mvij(p{x) = V{\xi — Thus 

j=l i<j 
RniP^^iO) ® Sn)R: = Di/„(0) . (37) 

Now we define 

for all t > 0, and we use the same notation to denote the extension of Z)i/„(t) to a density 
operator defined on all of Fe- From (jHTj) it may be seen that Di/„(t) is the solution of the von 
Neumann equation 



d 

dt' 

-Di/„(0) = Rn{P^»n (g) Sn)R* 



Z-Di/„(t) = [Hy^,Dy^it)] 



where Hi/n is the Hamiltonian ()23j) on Fe with L = — |A (E) /f2(j^) and V denoting the multi- 
plication operator 

{V(j)){xJ,y,k) = V{\x -y\)(t){x,j,y,k) 

on (L2(R3) ® £2(^N))®^ — a slight abuse of notation. 
Set X = 1/n. According to Statement 2 of Theorem 16.21 



A — *-0 



for all t > 0, where Afxi(t) = Afi{Di/n(t)) and satisfies the TDHF equation (j^ . One may 
verify that the single-particle density operator A/i(-Di/„(t)) equals Pn(t) ® Pn where P„ denotes 
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the orthogonal projector onto the span of {ci, . . . , e„} C i'^(N) and Pnit) is as defined in (j27j). 
On the other hand, we claim that = C,x(t) ® P„ where ^a(^) denotes the operator on 

L^(]R^) whose integral kernel satisfies 

d 



+ / {Vi\x-z\)-Vi\y-z\))Uz,z,t)dz Ux,y,t) 

- - / {Vi\x - z\) -V{\y - z\))Ux,z,t)U^,y,t)dz 
^xix,y,0) = 'il){y)i){x) = p„(x,?/, 0). (38) 
Accepting this claim for now, and noting that ||A/i(-Di/„(t))||i is always equal to n, we find that 

lim ||p„(t)-^;,(t)|| = (39) 

n^oo " " 

for_alH > 0, since ||A/'Ai(t) -Fi/„(t)||^ = \\{pn{t) - ^x{t)) ^ Pn\\^ = n \\pn{t) - ix{t)\\^- Equation 
is a small perturbation of (j^Hj) when n is large, and we may verify that ^x{t) converges to 



p*{t) by applying Gronwall's inequality to 

^^(P* - i\) = [L, p* - 6] + [V, (P* - 6) ® P*]:l 

+ [V, 6®(P*-6)]:l + llV, (6®G)f/]:l • (40) 

First, we get rid of the term [L, — ^x] by passing to the "interaction picture" and rewriting 
(jin|l as an equation for exp{itL){p^: — ^a) exp(— itL). Upon integrating and taking the trace 
norm one obtains 

\\p*{t) - Ut)\\, < M\V\\ f \\p.{s) - Us)\\,ds + -\\V\\ . 

Jo 

Gronwall's inequality implies that ||^A(i) — P*(i)||]^ = 0{l/n) as n — > oo for fixed t, which 
implies with that ||pn(^) — P*(i)||j^ — ^ 0, as asserted in the corollary. 
Finally, we verify that ^a(^) ® Pn satisfies when ^a(^) satisfies (jHH|) : 

= [L,F^/^] + iTrace3,4[^, {^x^Pnf'] 

- -Trace3,4 [V, (G ® Pnf'U^um)] (41) 

Th 
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where Trace3^4 denotes the partial trace over the third and fourth factors of the tensor product, 
and f/(i3)(24) denotes the unitary operator on ® that exchanges the first and 

third factors simple tensor products as it exchanges their second and fourth factors. The second 
term on the right-hand side of ()41|) is multiplied by l/ra to compensate for the Trace(P„) = n 
due to the partial trace over the fourth factor. But (j^T|) is the differential form of (|^^ for 
A = l/n. □ 
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